AlAs 2D electrons in an antidot lattice: Electron pinball with elliptical Fermi contours 
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We report ballistic transport measurements in a two-dimensional electron system confined to an 
AlAs quantum well and patterned with square antidot lattices of period a =0.6, 0.8, 1.0 and 1.5 /im. 
In this system two in-plane conduction-band valleys with elliptical Fermi contours are occupied. 
The low-field magneto-resistance traces exhibit peaks corresponding to the commensurability of 
the cyclotron orbits and the antidot lattice. From the dependence of the position of the peak 
associated with the smallest commensurate orbit on electron density and a, we deduce the ratio of 
the longitudinal and transverse effective masses mi/mt — 5.2 ± 0.4, a fundamental parameter for 
the anisotropic conduction bands in AlAs. 

PACS numbers: 72.20.-1, 73. 23. Ad, 75.47.Jn 
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The crystal structure and lattice constant of AlAs are 
closely matched to that of GaAs and, as a result, AlAs is 
widely used as a barrier material in AlAs/GaAs hetero- 
structures. However, AlAs differs greatly from GaAs 
in its band-structure. In contrast to GaAs where elec- 
trons occupy a single, isotropic, conduction-band min- 
imum at the center of the Brillouin zone, electrons in 
AlAs populate anisotropic conduction-band minima (val- 
leys) centered at the X-points of the Brillouin zone. The 
constant energy surface in (bulk) AlAs consists of six, 
anisotropic, half-ellipsoids (three full-ellipsoids). It is 
also possible to confine electrons to a modulation-doped 
AlAs quantum well, thereby forming a high-mobility two- 
dimensional electron system (2DES)fi If the quantum 
well is wider than about 5 nm, the electrons occupy 
only the two conduction-band valleys with their ma- 
jor axes in the plane, along the [100] and [010] crystal 
axes. We refer to these two valleys as X and Y, re- 
spectively [Fig. [TJb)]. Here we report ballistic transport 
measurements in a square anti-dot (AD) array^""^'^ pat- 
terned in a high-mobility AlAs 2DES. The data provide 
a direct measure of the anisotropy of the Fermi contour, 
or equivalently, the ratio of the longitudinal and trans- 
verse electron effective masses, a fundamental parameter 
of the AlAs conduction-band. Recently it has been shown 
that this anisotropy can be exploited to realize a sim- 
ple " valley- filter" device using a quantum point contact 
structure,^ Such a device may play an important role in 
"valleytronics" or valley-based electronics applications?^ 
or for quantum computation where the valley state of an 
electron might be utilized as a qubitf^ 

We performed experiments on 2DESs confined to a 
high-quality, modulation-doped, 15 nm-wide, AlAs quan- 
tum well grown by molecular beam epitaxy on a (001) 
GaAs substratCfi The AlAs quantum well is flanked by 
undoped layers of Alo.4Gao.6As on each side, and the 
doping is via a delta-layer of Si on the surface side at a 
distance of 100 nm from the well and 32 nm from the 
surface. We patterned a Hall bar sample, with the cur- 
rent direction along [100], using standard optical pho- 
tolithography. We then deposited a layer of polymethyl 
methacrylate (PMMA) and patterned the AD arrays us- 
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FIG. 1: (a) Micrograph of the AD lattice region with period 
a=0.8 fj,m. (b) The Fermi contours of AlAs in-plane valleys 
X and Y in fc-space. The Fermi wave vectors kF,x and kF,Y 
are indicated, (c) The first four commensurate orbits for the 
X and Y valleys that give rise to peaks in magneto-resistance 
for current along the a;-direction; these orbits have diameters 
in the y-direction that are equal to a multiple integer of the 
AD lattice period (see text). 



ing electron beam lithography. The PMMA layer served 
as a resist for a subsequent dry etching process used to 
define the AD holes. We used an electron cyclotron reso- 
nance etching system with an Ar/Cl2 plasma?^ at an etch 
rate of ~55 nm/min, to obtain small feature sizes without 
a degradation of the 2DES quality. The AD pattern was 
etched to a depth of ~80 nm, thus stripping the dopant 
layer and depleting the electrons in the AD regions. The 
micrograph of a section of one of our AD arrays is shown 
in Fig. [IJa). Each AD array is a square lattice and cov- 
ers a 20 /im X 30 fim area. There are four regions of AD 
lattice with different periods: a = 0.6, 0.8, 1.0, and 1.5 
/im as schematically shown in the inset of Fig.[2l The as- 
pect ratio d/a oi each AD cell is ~1:3, where d is the AD 
diameter. Finally, we deposited a front gate, covering 
the entire surface of the active regions of the sample to 
control the 2DES density. Following an initial back-gate 
biasing and brief illumination,^ we used the front gate to 
tune the total density (n^) from 2 to 5xl0^^/cm^. The 
density was determined from both Shubnikov-de Haas 
oscillations and Hall coefficient measurements that agree 
with each other. From the measurements on an unpat- 
terned Hall bar region in a different sample but from the 
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FIG. 2: Low-field MR traces for all four AD regions with 
periods (from top to bottom) equal to a = 0.6, 0.8, 1.0, and 
1.5 /xm. Traces are offset vertically for clarity, except for the 
bottom one. The Hall bar with the different AD regions is 
schematically shown on the right. 



FIG. 3; Low-field MR traces for the AD region with period 
a = 0.8 /im for Vg = -0.1 V to 0.15 V (from top to bottom), 
corresponding to a linear variation of the density ut from 2.27 



to 3.53x10' 
bottom one. 



/cm . Traces are offset for clarity, except for the 



same wafer, we obtain a mobility of ^10 m^/Vs at a typ- 
ical density of 3xl0"/cni2 and T=0.3 K. This gives a 
typical mean-free-path of '--^l /um. 

Figures [2] and [3] summarize our main experimental re- 
sults. Figure [2] shows the low-field magneto-resistance 
(MR) traces, measured as a function of perpendicular 
magnetic field (B), for all the AD regions. We observe 
two peaks, A and B, which are symmetric with respect 
to 5=0. Peak A, whose position is higher in field than 
peak B, is seen in the traces from all the AD regions. 
Peak B, on the other hand, is not observed in the a=1.5 
/im trace. In general, we observe that, as the AD lattice 
period becomes smaller, the positions of both peaks A 
and B shift to higher field values (as indicated by the 
dashed Hues). Figure [3] captures the gate-voltage (Vg) 
dependence of the MR traces for the a=0.8 /J,m AD re- 
gion. As we increase Vq to increase the 2DES density, 
peak A shifts to higher field values while peak B does 
not appear to shift. We have made similar observations 
in the other AD regions as Vg is varied. 

In order to analyze the data of Figs. [5] and [31 we first 
briefly review what is known about ballistic transport 
in AD arrays for GaAs 2DESs where the Fermi contour 
is isotropic. Low-field MR traces for such systems typi- 
cally exhibit commensurability peaks at magnetic fields 
where the classical cyclotron orbit fits around a group of 
ADs.-»24 Although there are subtleties associated with 
the exact shape of the AD potential and also the possi- 
bility of chaotic orbits that bounce from one AD bound- 
ary to another, the peak observed at the highest mag- 
netic field corresponds to the shortest period that fits 
around the smallest number of ADs. For an isotropic 
Fermi contour, this corresponds to a circular orbit, with 
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FIG. 4: MR obtained from numerical simulations (smooth 
curves are guides to the eye). Vertical lines indicate the ex- 
pected positions of the peaks for orbits XI, X2, X3 and Yl, 
based on Fig. [IJc) . We assumed equal densities for the two 
valleys and a current along the a;-direction. Bo is the magnetic 
field of the first commensurate orbit if the Fermi contour were 
circular. Inset: Simulation snapshot showing various types of 
trajectories: chaotic, pinned and skipping orbits, for X-valley 
electrons when XI peak occurs {B — mi/mt Bo). 



a diameter equal to the AD period, encircling a single 
AD. There have also been studies of ballistic transport 
in 2DESs with isotropic orbits in an anisotropic (rect- 
angular) AD latticeJ^ Experimental resultsj^i followed 
by theoretical analysisji^ indicated that the commensu- 
rability peaks are observed only when the orbit diam- 
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eter matches an integer multiple of AD lattice period 
along the direction perpendicular to the currentJ^ Based 
on the above considerations, we can predict the first four 
(smallest) commensurate orbits of the X and Y valleys 
that may give rise to MR peaks in our system; these are 
shown in Fig. [TJc). 

To evince this conjecture, we performed a kinematic, 
numerical simulation for our system, a 2DES with el- 
liptical Fermi contours in an isotropic AD lattice in the 
presence of a perpendicular magnetic field B. The sim- 
ulation details are similar to those described in Ref. [T^ . 
The AD boundaries are represented as hard-wall poten- 
tials where the electrons are scattered elastically upon 
collision. A set of typical electron trajectories are shown 
in the inset of Fig. H) Using the Kubo iormulsM^ we 
calculated the conductivity tensor and then obtained the 
resistance through its inversion. We performed the sim- 
ulations for two separate cases, one for X and one for 
y- valley electrons and present the results in Fig. U) The 
resulting MR traces indeed show peaks at or near the ex- 
pected values, namely orbits XI, X2, X3 for the X-valley 
and Yl for the F- valley [Fig. [^c)]. 

Now we focus our attention back on the experimental 
results. We associate peak A in our data of Figs. [2] and 
[3] with the shortest orbit XI in Fig. \V[^c) . From the field 
position of this peak, and if we assume that the electron 
density of the X valley is half the total density, we can 
directly obtain a value for the anisotropy (ratio of the ma- 
jor to minor axes diameters) of the elliptical orbits in our 
system, thus obtaining the effective mass ratio mi/mt^ 
However, there is a finite imbalance between the X and 
Y valley densities in our sample. Such imbalances can 
occur because of anisotropic strain in the plane of the 
sample and are very often present in AlAs 2DESs.^-^^ 
Therefore, we present here an analysis to determine the 
mi/mt ratio independent of the density imbalance. 

Consider a primary, commensurate orbit whose di- 
ameter in the direction perpendicular to the current is 
equal to the AD lattice period [orbits XI and Yl in 
Fig. [Ijc)]. These would give rise to MR peaks at fields 
Bp — 2hkp/ea where fc^ is the Fermi wavevector along 
the current direction. For the X and Y valleys, these 
wavevectors are kp^x and kp^Yi respectively, as shown in 
Fig. [IJb). For an elliptical Fermi contour, they are re- 
lated to the densities of the X and Y valleys, nx and 
ny, via the following relations: 

^'f,x = 27rnx \/mi/mt, kpy = '^'^ny \/mt/mi . (1) 

Note that the total density nT=nx+nY and the valley 
imbalance An—nx—ny- We can obtain nx from the 
Shubnikov-de Haas oscillations of the MR at high mag- 
netic fields or from a measurement of the Hall coefficient. 
Now consider orbit XI as shown in Fig. ^c) . Its associ- 
ated MR peak position Bp xi is given as: 

B^,x, = ^.[^{nT + An). (2) 

We use this expression to analyze our data. 
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FIG. 5: Summary of the density dependence of Bp for all 
four AD regions; Bp^a is the position of peak A observed 
in MR traces. The straight lines are linear fits using Eq. [2] 
The error bars reflect the uncertainty in the peak positions 
which were determined by subtracting second-order polyno- 
mial backgrounds from the MR traces. Inset: Slope (/?) of the 
Bp ^ vs ut lines of the main figure are plotted as a function 
of a~^. The dashed-line is a linear fit to the data. 

We assign peak A in our data to orbit XI and plot the 
square of its field position -Bp^ as a function of nx in 
Fig. [51 It is clear that for all four AD lattice regions, Bp ^ 
varies linearly with np as expected from Eq.[21 Moreover, 
we obtain the slopes {(3) and the intercepts of the lines in 
Fig. [5] by performing a least-squares fit of each data set. 
Note that according to Eq. [11/3 = h-^y^mi/mt/ire'^a'^, 
and the intercept is equal to /3AnJ^ Finally, we plot /3 
as a function of in the inset of Fig. \5\ This figure 
shows that, consistent with the prediction of Eq. [H f3 
indeed depends linearly on and the line has a zero 
intercept. 

From the slope, A/3/A(a^^), of the line in Fig.[5linset, 
we can deduce the effective mass anisotropy ratio: 

mi/mt = TT^e^h^ [Al3/A{a-^)] ' . (3) 

Note that this mass anisotropy ratio is related to the 
slope of the line in Fig. [5l inset by a pre- factor containing 
only physical constants {n'^ e* / h'^) . Our data analysis and 
determination of this ratio are therefore insensitive to pa- 
rameters such as density imbalance between the two val- 
leys. From data of Fig. [51 we obtain mi/mt = 5.2 ±0.4, in 
very good agreement with the ratio mi/mt = 5.2±0.5 de- 
termined from recent ballistic transport measurements in 
AlAs 2DESs subjected to one-dimensional, periodic po- 
tential modulations.l2 Given the cyclotron mass in AlAs, 
m*—^mimt — 0.46 mor^ we deduce mi — (1.05±0.1)too 
and nit = (0.20±0.02)mo, where niQ is the free electron 
mass, in good agreement with the results of the major- 
ity of theoretical and experimental determinations of the 
effective mass in AlAsJ^ 
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FIG. 6: Summary of the density dependence of Bp of peak A 
and B for AD region: (a) a = 0.6 /im and (b) a = 0.8 /im. The 
lines marked XI are linear fits to peak A using Eq. (2] Lines 
X2, X3 and Yl are the predicted peak positions calculated 
using equations similar to Eq. [2](see text). 

Several other features of the data presented here are 
noteworthy. From the intercepts of the hnear fits in Fig. [5] 
we can determine the valley density imbalance for each 
AD lattice. Such analysis gives An — 1.5, 1.2, 1.0, and 
1.3 xlO" /cm2 (±0.2x10" /cnP) for the AD regions 
with a — 0.6, 0.8, 1.0, and 1.5 /im, respectively. Such 
a variation of valley imbalance for different AD regions 
may come from non-uniform residual strain across the 
sample. Note that, because of the close proximity of the 
different AD lattice regions, we expect this variation to 
be small, consistent with the An values deduced from the 
above analysis. 

As for peak B, it is tempting to associate it with orbits 



X2 or Yl in Fig. [TJc). This is qualitatively consistent 
with the data of Fig. [5] which indicate that peak B moves 
to smaller values of magnetic field as the period of the AD 
lattice is made larger. Moreover, peak B becomes weaker 
with increasing AD lattice period and disappears for the 
largest period a = 1.5 /im. This is also consistent with 
the larger size of the X2 and Yl orbits (compared to the 
XI orbit), and the fact that for a = 1.5 /im, the lengths of 
these orbits become large compared to the electron mean- 
free-path. Quantitatively, using the values oimi/mt and 
An obtained above, we can modify Eq. [^and determine 
the expected peak positions associated with the X2, X3 
and Yl orbits »2£ As illustrated in Fig.[6l we find that the 
predicted peaks for X2, X3 and Yl orbits are close to each 
other in field and approximately straddle the observed 
positions of peak B in Fig. [21^^ It is possible then that 
peak B may originate from a superposition of X2, X3 and 
Yl peaks that cannot be resolved in our experiment. We 
cannot rule out, however, that peak B may be strongly 
influenced by non-linear orbit resonances in the system. 
Such resonances are known to occur for orbits with long 
trajectories in the presence of a smooth AD potential^i 

In conclusion, we performed ballistic transport exper- 
iments in AD lattices imposed on an AlAs 2DES where 
the electrons occupy two valleys with anisotropic Fermi 
contours. The low-field MR traces exhibit two sets of 
peaks. From the analysis of the positions of the peak 
associated with a commensurate orbit with the shortest 
trajectory [XI orbit in Fig.[ljc)], we deduced the effective 
mass anisotropy ratio mi/mt, a fundamental parameter 
of the AlAs conduction-band dispersion that cannot be 
directly measured from other transport experiments. 
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To calculate the peak positions for orbit X2 and X3, simply 
replace a with 2a and 3a in Eq. [2] For orbit Yl, replace 

mi /mt with mt/mi and An with — An. 
In Fig.[S]the expected Bp for Yl orbits are lower than Bp 
for X3 because ny < nx in our experiments. 



